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Abstract 

We propose an interesting mechanism to reconcile the recent experiments of the 
Michelson-Morley type and slowdown of the velocity of light in dark matter with a 
fractional electric charge when the index of refraction of dark matter depends on the 
frequency of a photon. After deriving the formula for the velocity of light in a medium 
with the index of refraction n{uj) in a relativistic regime, it is shown that the local 
anisotropy of the light speed is proportional to the second order in n(uj) — 1. This result 
implies that the experiments of the Michelson-Morley type do not give rise to a stringent 
constraint on the slowdown of the velocity of light in dark matter with electric charge. 
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1 Introduction 



It is nowadays well known that a large portion, about 23 percentages, of the total energy 
density in the present universe is occupied by unknown material, what we call "dark matter". 
The particles, of which dark matter is made, are usually thought to be weakly interacting 
massive ones, in other words, WIMPs, which are absent in the Standard Model (SM). These 
hypotheoretical particles are supposed to have mass of order ranging from GeV to TeV, 
and cross section of their annihilation process is comparable to that of the weak interaction. 
Thus they do not completely annihilate in the course of evolution of the universe and their 
mass density in the present universe could be of the order of the critical density pc ~ 0.52 x 
l{i~^GeV/cm^ (more precisely, pou ~ 1-2 x lQ~^GeV / cm^ averaged over large distance scales) 

[!]• 

The difficulty of terrestrial experiments for detecting dark matter by direct or indirect 
methods makes it impossible to clarify its physical properties and therefore answer various 
questions associated with dark matter such as the mechanism of the generation in the early 
universe. On the other hand, the density of dark matter in clusters of galaxies is determined 
by various methods of measurement of the gravitational potential. For instance, mass distri- 
bution in a cluster is obtained by the method of gravitational lensing. The result is that most 
of the mass is due to dark matter distributed smoothly over the cluster. Put differently, dark 
matter, like the usual matter, is more dense in galaxies and it appears as if luminous matter 
is embedded into the cloud of dark matter of larger size-galactic halo. It is easy to understand 
this concentration of mass of dark matter near clusters of galaxies from a physical viewpoint 
since the gravitational interaction attracts dark matter near the clusters compared to empty 
outer space. 

Given our current ignorance of the dark sector in the universe, it seems prudent not to re- 
strict our interests only in WIMPs, though they must be most plausible, and open our mind 
to many possibilities. Indeed, though the particles consisting of WIMPs carry no electric 
charge, the notion of charged massive particles (CHAMPs) has also appeared periodically 
[2]- [6]. Since CHAMPs, which normally carry a full unit of electric charge, receive stringent 
constraints from searches for heavy hydrogen to direct detection in the underground experi- 
ments, the possibility of CHAMPs carrying a fractional charge has been recently investigated 
[7, 8]. 

In this article, we shall assume that the earth itself is embedded into a medium of dark 
matter with a fractional electric charge. Then, it is quite natural to think that the existence 
of such a dark matter manifests itself the index of refraction n{u) depending on the (angular) 
frequency of light since the dark matter can interact with the gauge field in a direct way 
[9].^ Moreover, we can ask ourselves if the recent experiments of the Michelson-Morley type 
[11, 12] shed some hght on the possibility of detecting the dark matter on the earth through 

^Precisely speaking, an electrically- neutral particle propagating through matter is also refracted if it couples 
to a charged particle even indirectly. For instance, a neutrino can couple to an electron via weak interaction 
and its refractive property has an important effect on neutrino propagation when the neutrino is massive [10]. 



1 



measuring the local anisotropy of velocity of light or not. 

One of motivations behind the present study comes from a recent claim [13] that in case 
of the OPERA superluminal neutrinos [14], such an assumption shows a clear tension and 
might be in conflict with the most recent experiments of the Michelson-Morley type about the 
local anisotropy of velocity of hght [11, 12]. In order to relax this tension, it is useful to recall 
that the index of refraction of a medium is in general dependent on the frequency of light. 
If we take account of this fact together with its property under the Lorentz transformation 
in a proper manner, a term depending on the frequency provides an additional contribution 
to the velocity of light and then the resulting constraint would become consistent with the 
experimental results of the Michelson-Morley type. 

This article is organized as follows: In the next section, we derive a formula of propagation 
of light in moving media based on special relativity. In Section 3, it is shown that using this 
formula the local anisotropy of the light speed is proportional to the second order in n{uj) — 1. 
Section 4 is devoted to discussion. 



2 Velocity of propagation of light in moving media 

In this section, we review special relativity [15], in particular, on the basis of both the addition 
law of relativistic velocities and the invariance of phase of a plane wave, we shall derive the 
velocity formula for propagation of light in moving media. 

It is known that the Lorentz transformation for inertial coordinates can be used to derive 
the addition law of relativistic velocities. To show this fact explicitly, let us consider a 
physical setup where we have a moving point P whose three-dimensional velocity vector u' 
has spherical coordinates {u', 6', 0') in the inertial frame K', and the frame K' is moving with 
velocity v = q/3 in the direction of the xi axis with respect to the inertial reference frame 
K. Here q is a universal limiting speed, in other words, the maximum speed of all physical 
entities [16] -[21]. 

Then, one wishes to calculate the velocity u of the point P as seen from the inertial 
frame K. To do so, let us first note that the Lorentz transformation takes the form for the 
differential expressions 

dxo — 'y{dxQ + f3dx[), 
dxi — ^{dx'i + Pdx'o), 

dx2 — dx'2, dxs — dx'^, (1) 
where we have defined 7 = , ^ ^ — . Since the components of velocity are defined by u[ ~ 

and Ui — ci^, they transform as 

Ui = r, U2 = —, Us — —. (2) 
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For later convenience, let us rewrite this transformation to a more general form 

where u\\ and tix indicate the components of velocity parallel and perpendicular to v, re- 
spectively. Moreover, it turns out that the magnitude u oi u and its polar angles in the two 
inertial frames are related as follows: 

(j) = 0', 

^ u' sin 9' 

tanO - 



7(m' cos 6*' + f)' 

Ju'^ + v^ + 2u'v cos 6' - (H^^saflji 

u = ^ ^' ■ (4) 

1 + ^cos^' ^ ^ 

The inverse transformation for u' in terms of u can be easily obtained by interchanging primed 
and unprimcd quantities and simultaneously changing the sign of v. 

In what follows, for simplicity, let us focus on the case of the parallel velocities, 9' = 0. In 
this case, the magnitude of u takes the form 

~^ 

If one particularly sets u' = ci, this expression reads u = q, which simply means the postulate 
of a universal limiting speed, which is an alternative fundamental principle for the principle 
of invariant speed of light in special relativity [16, 17, 18, 19]. 
Next, let us choose = thereby giving rise to 

n(aj') 

^11 = 1 ^ ' (6) 

^ n(oj')ci 

where uu' is the (angular) frequency of light in the inertial frame K' and n{u!') is the index of 
refraction of a medium. The important point is that the index of refraction of a medium in 
general depends on the magnitude of the frequency of the photon. 

In order to determine n{uj'), we make use of the fact that the phase of a plane wave is 
an invariant quantity under the Lorcntz transformation since the phase can be identified 
with the mere counting of wave crests in a wave train, an operation that must be the same 
in every inertial frame 

ip = kf.x'^ = ujt-k-x = u't' - k' ■ x'. (7) 

With the frequencies a; = q/cq and a;' = q/cq, the Lorentz transformation for the wave-number 
vector k^ reads 

/cq = 7(A:o-/3-fc), 
Ml = 7(^11 -^^o), 

^1 = k^. (8) 
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For light waves, \k\ = ko and \k'\ = /cq, so using 9 = coming from 9' = via Eq. (4) we have 
the relation 



u' = 7(1 
lolds wht 

transforms to first order in ^ as 



1 + 



(9) 



where the last relation holds when B = - is small. Thus, the index of refraction of a medium 

Cl ' 



n(uj 



n{uj) — f5u} 



dn 



(10) 



Plugging this expression into Eq. (6) (and taking ±v into consideration at the same time), it 
is easy to show that for medium flow at a speed v parallel or antiparallel to the path of the 
light, the speed of the hght, as observed in the laboratory (X- frame), is given to first order 
in V by 



u 



± 



Cl 



n{u) 



±v 



1 



1 ^ uj dn{uj) 



n(uy n(uy duo 



(11) 



where superscripts + and — on M|| denote the speed parallel and antiparallel to the light 
path, respectively. Formula (11) is the main result of this section and will be used in the next 
section. 



3 The recent experiments of Michelson-Morley type 

We wish to apply the formula (11) to the problem that the recent experiments of the 
Michelson-Morley type could or could not give us useful information on detection of dark 
matter with electric charge. 

Before doing that, let us present some ideas behind our theory. Search for dark matter is 
underway, but no conclusive evidence has been obtained thus far. A direct way for the search 
is carried out in experiments trying to detect energy deposition in a detector caused by elastic 
scattering of dark matter off a nucleus. Besides the direct search, there are experiments to 
search for dark matter particles indirectly, which include the search for products of dark 
matter annihilation. One of the most promising indirect ways is to search for monocnergetic 
photons which are emitted in two-body annihilation processes X+X — )■ 7+7, X+X — )> Z°+7 
where X, X describe a dark matter particle and its antiparticle, respectively. By the crossing 
symmetry, these processes read X-I-7— )-X + 7,X-|-7^X-|- Z^. If the photon energy is 
small compared to the weak energy scale (i.e., the energy of the photon beam is below the 
threshold for production of a single Z°), only the former process is allowed to occur and it is 
nothing but the elastic Compton scattering. 

At this stage, let us recall that there is a well-known connection between the index of 
refraction and the Compton scattering amplitude where the standard formula takes the form 
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[22] 



, , 27rcfNF(uj) , , 

n o; - 1 + 2 1 12 

where N is the number of scattering centers per unit volume and F{uS) is the forward Compton 
scattering amphtude for the photon scattering off the medium which is a function of the 
(angular) frequency a; of a photon. Put differently, in this approach the dark matter with 
electric charge plays a role of a medium for the photons and provides a dispersive effect on 
light propagation. 

Since we assume in this article that dark matter carries a fractional electric charge asso- 
ciated with [/(I) EM symmetry, it must be a stable Dirac particle with spin | like protons. 
Then, the amplitude for the forward scattering of a photon from a polarization state e\ to 
the one 62 by the dark matter must be of the general form [22] 

F(a;) = i)) [Fi(a;)e; • ei + %F2{u)d ■ (e^ x ei)] ^i, (13) 

where ipi{i/jf) is the wave function of the initial (final) dark matter and a is the spin matrix 
of the dark matter. In general, both Fi{uj) and F2(cj) are complex and have dispersive and 
absorptive parts. If one averages over dark matter spins in the amplitude one is left only 
with Fi{uj). The amplitudes ^1(0;) and ^2(0;) are separable if one can do experiments with 
polarized dark matters: Fi{uj) corresponds to parallel and ^2(0;) to perpendicular linear 
polarization vectors of the initial and final photons, respectively. Since we are interested in 
the spin-averaged forward amplitude, we will focus on only Fi{uj) henceforth. (We therefore 
define F = F^.) 

With the assumption of causality and analyticity for the forward scattering amplitude 
F{uj), the once-subtracted dispersion relation is given by 

ReF(.) = «eF(0) + H ^'^fP^,, (14) 

where P denotes the principal value and ouq is the threshold for producing a single dark matter 
which is approximately taken to be Uq ^ Mxcf ~ 100Ge\^. Using the optical theorem, the 
imaginary part of the forward elastic scattering amplitude is related to the total cross section 
a{uj) by 

a{uj) = ^ImF{u). (15) 

Furthermore, it is known that the forward elastic scattering amplitude at a; = is real and 
given by the Thomson formula 
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where ee and Aix are respectively the charge and mass of the dark matter particle. Thus, 
one can rewrite ReFiuj) in Eq. (14) to the form 



Substituting this expression into the standard formula (12), we obtain 

Ren{u:) - 1 = -j-^ + — P / dw' ' . 18 

Since Imn{u) does not appear any more and n{uj) in the previous section is in fact equivalent 
to Ren{uj) in this section, for simplicity we rewrite Ren{ijj) as n(u;) from now on. 

Now let us consider the case of the low photon energy of cu <C a;o.^ It is then sufficient to 
limit ourselves to an expansion up to the constant order in uj"^ ^ 

5{uj) ^ n{u) -l^% + 52, (19) 

where 5q, 82 are defined as 

2nN{eef , _ ciN ^ ^ X^') 

lu!o 



, P duj'^^. (20) 



Note that observations of 7-ray burst (GRB) emission from radio energies to tens of GeV [24] 
require us that the energy dependence of the speed of the photon is so small that 5q must be 
very tiny. Moreover, the recent ICARUS result [25] implies that 82 are very small as well.^ 

However, this smallness is not a real problem for the present purpose. The key point is 
that if we regard dark matter as a kind of usual medium for light, the frequency dependence 
must be always given by (19) at least for the low energy photons. Then, an important relation, 
which will be utilized shortly, can be obtained 

a;^ = -2(n-l-5.). (21) 

Next, following the procedure in Ref. [13], let us calculate the local anisotropy of the 
speed of light whose experimental upper bound has been already obtained in the experiments 
of the Michelson-Morley type [11, 12] 

—\ ^Ix 10-1^. (22) 

C J Exp 



^In the most recent Michelson-Morley experiments [11, 12], photons with the order IQQTHz are utihzed, 

which corresponds to the order leV. 

^Expanding in powers of w is known to be a useful tool for obtaining various interesting low-energy cross 

sections in the analysis of the Compton scattering [23] . 

^Although E = huj corresponds to the energy of a photon, we will regard it as a typical energy scale 

controlling the system under consideration, in other words, the energy of the ICARUS neutrinos. 
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In these experiments, the parallel average velocity is measured 



T^. (23) 
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Using Eqs. (11) and (21), this quantity can be evaluated to be 

^11 = n - Q ■ ^24) 

Then, the theoretical local anisotropy of the speed of light is calculated to be 

^ = ito.4(^)%,^«4f^V.^ (25) 



ClJ \CiJ 



where we have used 5 a; 0,/3 <^ 1 and S !^ 82. It is worthwhile to notice that compared to 
the result in Ref. [13] which is proportional to the hnear order in 5, our result (25) becomes 
proportional to the second order in S since we have taken account of the frequency dependence 
and the Lorentz transformation of the index of refraction n{u). 

In order to show concretely that Eq. (25) yields a very weak constraint on the local 
anisotropy of the speed of light, let us apply it to the OPERA result [14] together with the 
SN1987A result [26]. Here note that the OPERA result has been recently refuted by the 
ICARUS group [25]. The reason why we use the OPERA result is two-fold. First, Ref. [13] 
has discussed the OPERA result and shown that there is some tension between the recent 
Michelson-Morley experiments and the models based on CHAMPs owing to a large value of 
S in the OPERA result. But we wish to show that even in this large value of 6 there appears 
no tension between them since (25) becomes proportional to the second order in 6 in our 
calculation. And the second reason is that even if the OPERA main result, in particular the 
size of S might be wrong, the other result, i.e., 5 is independent of the neutrino energy, is 
consistent with that of GRB emission [24] and seems to be true. 

The result of SN1987A implies that the velocity of a neutrino is almost the same as a 
universal limiting speed 

v^KiCi. (26) 
Then, we can obtain the following relation 



a r, - £L 



n _J' n 

n n 



2.37 X 10-^ (27) 



where we have used the OPERA results in the last step. The OPERA result suggests that 6 
is independent of the neutrino energy, so this equation implies 

62^6^ 2.37 X 10"^ (28) 
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Substituting this value and the equatorial rotation speed of the earth at the observed place, 



/3 = ^ = 1.55 X 10 ^ [13] into Eq. (25), one arrives at the result 

-] ^ 4.6 X 10~^^ (29) 



Cl 



c J 



OPERA 



Note that this value is much smaller than the experimental upper bound (22), so it is not 
in conflict with the recent experimental results of the Michelson-Morley type. Of course, 
our result docs not refute the result obtained in Ref. [13] directly since in our approach at 
hand the index of refraction depends on the frequency of a photon whereas in the approach 
considered in Ref. [13] the index of refraction is assumed to be independent of the photon 
frequency. Actually, when dark matter does not carry an electric charge like a neutrino, Eq. 
(20) means Sq = 0, thereby making it impossible to get Eq. (21) owing to = 0. 

For completeness, let us move on to a general case of the arbitrary photon energy. In this 
case, we can proceed the argument in a perfectly similar way to the case of the low energy 
photons. The relation (21) is now replaced with 

u^^-2{n-l-~5), (30) 
auj 

where S is defined as 

s^siUpr^.^^^, (31, 

Then, it turns out that local anisotropy of the speed of light is exactly evaluated to be 

C \CiJ v? 

With the reasonable assumptions as before, 5 = n — 1~0,/3^1 and (5 ~ we can obtain 
a similar result to Eq. (25) 

^^a('-\"~5^^4'-]\\ (33) 



C \CiJ \CiJ 

Compared to the low energy case, a slight modification here appears in the 5 which is depen- 
dent on the photon frequency as can be seen in (31). However, except this fact, the essential 
feature remains unchanged so the conclusion obtained in the case of the low enery photons 
still holds even in this general case . 



4 Discussion 

Since the introduction of dispersion relations into elemetary particle physics, originally within 
the context of quantum field theory [22], a large number of literatures have grown up on their 
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theoretical basis, on extensions and applications to new processes, and on their comparison 
with experiment with great success. The most advantageous point is that dispersion relations 
are very universal in the sense that they are formulated only on the basis of fundamental 
principles of quantum field theory such as causality, the Lorentz invariace and analyticity. 

As a simple application of the dispersion relations, in this article, we have studied the 
velocity of light in dark matter with a fractional electric charge and found that the recent 
experiments of the Michelson-Morley type [11, 12] do not provide a stringent condition on 
the slowdown of the velocity of light if the index of refraction of dark matter depends on the 
frequency of a photon. The experiments of the Michelson-Morley type have played a critical 
role in proving the non-existence of the aether, i.e., a medium that was once supposed to 
fill all space and to support the propagation of electromagnetic waves, and might reveal a 
breakdown of the Lorentz invariance at the Planck scale in future. But they do not seem to 
shed some light on the change of the velocity of light in dark matter with charge. Perhaps, 
if the ICARUS result is correct, this situation will not change even in the future experiments 
of the Michelson-Morley type.^ 

Nevertheless, the experiments using the photons will provide important insight into the 
properties of CHAMPs in future. For instance, polarized Compton scattering will yield in- 
formation on spin-structure of the dark matter. Moreover, if one increases the energy of the 
incoming photon beam, one can probe the internal structure of the dark matter [23]. 

Finally, we wish to close this article by mentioning two future problems to be solved. An 
interesting application of the refractive effect in our model is calculation of the force acting 
on a macroscopic body when dark matter flux passes through. Simple calculation turns 
out to give a force proportional to \n{(jj) — 1| which could be measured in the Eotvos-Dicke 
experiment. Furthermore, the refractive effect in a high-temperature background such as in 
the early universe deserves study in future. It is of interest to notice that the presence of the 
electric charge might ensure the stability of the dark matter from another angle. 
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